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ABSTRACT 


The  purpose  of  this  thesis  is  to  model  digital  communications  systems  in  the  time 
domain  using  MATLAB  Simulink  and  the  Communications  Toolbox  as  well  as  to 
determine  and  verify  system  performance  in  the  presence  of  additive  noise  and  co¬ 
channel  interference.  While  the  theoretical  results  are  available  for  the  effect  of  wideband 
gaussian  noise  on  the  performance  of  digital  communications  systems,  determining  the 
performance  of  a  system  in  the  presence  of  noise  and  co-channel  interference  is  best  done 
by  computer  simulation.  Time  domain  modeling  allows  the  visualization  of  the 
communication  signal  at  various  stages,  and  “Monte  Carlo”  type  simulations  establish  the 
bit  error  rates  under  realistic  conditions  of  noise  and  co-channel  interference  as  well  as 
for  different  transmitter/receiver  parameters,  different  channel  parameters,  and  different 
types  of  interfering  signals.  It  should  be  noted  that  the  simulation  does  not  account  for 
system  non-linearities. 
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1.  INTRODUCTION 


A  military  site  or  a  platform  typically  has  a  large  number  of  communication 
systems  operating  simultaneously  in  both  transmit  and  receive  modes.  The  interference 
of  one  system  with  another  depends  on  available  frequency  allocations,  the  physical 
separation  of  transmit  and  receive  antennas,  non-linearity  of  transmitter  output  stages,  and 
so  on.  The  objective  of  this  study  is  to  obtain  a  computer-based  prediction  of  noise  and 
co-channel  interference  effects  on  the  probability  of  bit  error  for  M-ary  Phase  Shift 
Keying  (MPSK)  digital  communication  systems.  To  accomplish  this,  a  model  was 
developed  using  SIMULINK  and  the  Communications  Toolbox.  The  performance  of  the 
model  has  been  verified  by  comparing  the  simulation  results  and  the  theoretical  results  for 
the  bit  error  probability  in  the  presence  of  additive  white  Gaussian  noise  (AWGN).  The 
model  has  been  subsequently  modified  to  include  co-channel  interference,  and  the  bit 
error  probability  of  MPSK  with  AWGN  and  co-channel  interference  was  obtained  via 
simiilation. 
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II.  SIMULINK  AND  THE  COMMUNICATIONS  TOOLBOX 


SIMULINK  is  a  program  for  modeling  linear  and  non-linear  dynamic  systems  in 
the  time-domain.  Models  in  SIMULINK  are  represented  in  block-diagram  form  and  can 
be  assembled  from  block  libraries  and  sub-libraries.  Furtiiermore,  SIMULINK  allows  the 
display  of  results  and  changes  of  certain  model  parameters  without  intermpting  the 
simulation.  Since  SIMULINK  is  built  upon  the  MATLAB  numeric  computation  system, 
it  offers  direct  access  to  the  MATLAB  workspace  and  MATLAB ’s  mathematical  and 
engineering  functions.  SIMULINK’ s  main  features  are: 

•  modeling  and  analysis  of  dynamic  systems,  including  linear,  nonlinear, 
continuous,  discrete,  and  hybrid 

•  flexible  "open  system"  environment  that  allows  addition  of  new  blocks  to 
Simulink 

•  seamless  interface  with  MATLAB ’s  built-in  math  functions,  2-D  and  3-D 
graphics,  and  add-on  toolboxes  for  specialized  applications 

•  choice  of  methods  of  running  a  simulation  (menu-driven  on-screen  or  batch¬ 
mode) 

•  an  optimized  computer  platform  implementation  that  ensures  fast  and  accurate 
results 

•  unlimited  model  size. 

A  typical  SIMULINK  session  starts  by  either  defining  a  new  model  or  recalling  a 
previously  defined  model  and  then  proceeds  to  a  simulation  of  the  performance  of  that 
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model.  In  practice  these  two  steps  are  often  performed  iteratively  as  the  model  designer 
creates  and  modifies  a  model  to  achieve  the  desired  behavior  [Ref.l]. 

The  Communications  Toolbox  is  a  collection  of  MATLAB  functions  and 
SIMULINK  blocks  for  model  development  and  simulation  in  the  communications  area. 
The  functions/blocks  are  organized  in  the  following  sub-categories:  Data  Source,  Source 
Coding  and  Decoding,  Error-Control  Coding,  Modulation  and  Demodulation, 
Transmission  and  Reception  filters,  Transmitting  Channel,  Multiple  Access,  and 
Synchronization  and  Utilities  [Ref.l]. 
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m.  PHASE  SHIFT  KEYING 


Modulation  of  a  sinusoidal  carrier  allows  the  conveyance  of  information  in  binary 
form.  The  baseband  information  bit  stream  modulates  the  carrier  by  means  of  discrete 
changes  in  the  carrier  amplitude  (amplitude-shift  keying  or  ASK)  carrier  frequency 
(frequency-shift  keying  or  FSK)  or  carrier  phase  (phase-shift  keying  or  PSK)  or  some 
combination  of  the  above.  Each  different  state  of  the  carrier,  known  as  a  symbol, 
corresponds  to  one  or  more  bits  of  baseband  information  [Ref  2.]. 


A.  BINARY  PHASE-SHIFT  KEYING 


In  binary  phase-shift  keying  (BPSK),  the  phase  of  a  constant  amplitude  carrier 
signal  is  switched  between  two  values  according  to  the  two  possible  “messages”  mi  and 
m2  corresponding  to  binary  1  and  0,  respectively.  Normally,  the  two  phases  are  separated 
by  180  degrees.  If  the  sinusoidal  carrier  has  amplitude  Ac,  then  the  average  energy  per  bit 


=^A^Ti,,  and  the  transmitted  BPSK  signal  can  be  represented  as: 


^  BPSK  (f) 


^BPSK  W 


■if- 

■Jf' 

=  -(^ 

i  n 


0<t<T^ 

( binary  1) 

(1) 

0<t<Ti, 

( binary  0) 

(2) 

0<t<Ti, 

( binary  0) 

(3) 
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where  is  the  average  energy  per  bit,  is  the  bit  duration,  and  a  rectangular  pulse 


shape  p(t)=;i[(/  is  assumed.  For  this  signal  set,  the  single  waveform 


0<t<T, 

can  be  defined.  Using  this  signal,  the  BPSK  signal  set  can  be  represented  as 

^BPSK  ~  (f)  } 


(4) 


(5) 


B.  SIGNAL  CONSTELLATION  FOR  BPSK 


Figure  1.  Signal  Constellation  for  BPSK 

A  BPSK  signal  can  be  geometrically  represented  in  “  signal  space”  as  shown 
below.  Such  a  representation  is  referred  to  as  the  symbol  constellation  and  provides  a 
graphical  representation  of  the  complex  envelope  of  a  BPSK  signal  for  the  two  possible 
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symbols.  The  x-axis  of  the  constellation  diagram  represents  the  “in-phase”  component  of 
the  complex  envelope,  and  the  y-axis  represents  the  “quadrature”  component  of  the 
complex  envelope.  The  distance  between  the  signals  on  a  constellation  diagram  indicates 
how  different  the  modulation  waveforms  are  and  determines  how  well  a  receiver  can 
differentiate  between  all  possible  symbols  when  random  noise  is  present.  The  larger  the 
signal  distance,  the  better  the  chance  of  correct  symbol  detection. 


C.  PROBABILITY  OF  BIT  ERROR  FOR  BPSK 


An  important  measure  of  performance  used  for  digital  modulation  is  the 
probability  of  symbol  error,  .  It  is  often  convenient  to  specify  system  performance  by 
the  probability  of  bit  error  P^ ,  even  when  decisions  are  made  on  the  basis  of  symbols 
rather  than  bits.  The  relationship  between  Pg  and  Pg  for  orthogonal  signaling  is: 


(6) 


logjM 

For  the  BPSK  modulation  (M=2),  the  symbol  error  probability  is  equal  to  the  bit 
error  probability.  When  the  signals  are  assumed  equally  likely  and  signal  5;(r)(i=l,2)  is 


transmitted,  the  received  signal,  r{f),  is  equal  to  sXt)  +  n{f),  where  n(t)  is  modeled  as 
additive  white  Gaussian  noise  (AWGN).  The  antipodal  signals  (signals  of  equal 
amplitude  and  opposite  polarity)  Sy{t)  and  s^it)  are  (cf.  Eq.5): 

Sr(t)  =  ^MO  0<t<T,  (7) 

=  0<t<T,  (8) 
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The  decision  stage  of  the  detector  will  choose  the  5,-  (t)  with  the  largest  correlator  output 
Zi  (t) ,  or  in  this  case  of  equal-energy  antipodal  signals,  the  detector,  using  the  decision 
rule,  decides: 

5i(0  if  z(r)>)o 

S2  (t)  otherwise 

where  yo  denotes  the  decision  threshold  (equal  to  0  for  equally  probable  antipodal  signals) 
and  z  (7)  is  the  correlator  output  at  time  T.  Two  types  of  detection  error  can  be  made. 
The  first  type  of  error  takes  place  if  (t)  is  transmitted  but  the  noise  is  such  that  the 
detector  measures  a  negative  value  for  z(r)and  decides  (incorrectly)  that  signal  53(0  was 
sent.  The  second  type  of  error  takes  place  if  signal  ^2(0  is  transmitted  but  the  detector 
measures  a  positive  value  for  z(T)  and  decides  (again  incorrectly)  that  signal  s^(t)  was 
sent.  Therefore,  the  probability  of  bit  error,  is  the  sum  of  two  conditional 
probabilities: 

=  P(S2 1  S,)P(S,)  +  P(s,  I  S2)P(S2)  (9) 

For  the  case  when  the  symbol  probabilities  are  known  a  priori  and  the  symbols  are 
equally  probable  (which  is  mostly  the  case): 

P(;y,)  =  P(s2)=|  (10) 

The  expression  for  the  bit  error  probability  becomes: 

P,=^P(s2\s,)+^Pis,\s2)  (11) 
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Because  of  the  symmetry  of  the  probability  density  functions  (pdf  s)  of  the  sum  of  signal 
and  noise,  P(z  1 5i)  and  P(z  1 52)the  two  conditional  probabilities  are  also  identical  and 

P,,  =  P(52l5l)  =  P(5j^2)  (12) 

The  probability  of  a  BPSK  bit  error  Pg  is  numerically  equal  to  the  area  under  the  “tail”  of 
either  pdf  p{z  I  ^'i)  or p{z  I  s^)  that  falls  on  the  incorrect  side  of  the  threshold.  We  can 
therefore  compute  by  either  integrating  p(zi5i)  between  the  limits  -<»  and  or  by 
integrating  p{z\s2)  between  the  limits  jg  «>.  Hence, 


Pb=  °\piz\s2)dz  (13) 

where  the  conditional  pdf  sp(z  I  Sj)  (/ =  1,2)  are  Gaussian  with  mean  value  a,  , 
variance  Oo,  and  the  optimum  threshold,  jg,  is  The  probability  of  bit  error 

for  BPSK  is: 


J 

ro=(‘»i+a2)/2 


1  r  1  y-  2  ^2  \2 1 


exp[-^(^^)^]& 

2  (To 


(14) 


If  we  introduce  u=— — —  and  Of^U—dz, 

<5’o 

the  integral  simplifies  to 


Pb  = 


K=oo 

/ 


^  exp(r^)du  = 


u={ai-(ii)l2aQ 


2(T„ 


(15) 


where  <3g  is  the  standard  deviation  of  the  noise  at  the  output  of  the  correlator.  The 
function  Q{x)  is  defined  as: 


1  2 

Q{x)  =-f^\exp(-^)du 
dlTC  i  2 


(16) 
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For  equal  energy  antipodal  signaling,  such  as  the  BPSK  format  in  (5),  the  receiver  output 
signal  components  are  when  s^it)  is  sent  and  fl,  when  jjCOis  sent. 

For  AWGN  the  noise  variance  at  the  correlator  output  is  equal  to  Nq/2  so  that  we 
can  rewrite  probabihty  of  bit  error  as  follows  [Ref.3]: 


i— exp(-— 


(17) 


(18) 


D.  QUADRATURE  PHASE-SHIFT  KEYING 

In  quadrature  phase-shift  keying  (QPSK),  two  bits  are  transmitted  in  a  single 
modulation  symbol  The  phase  of  the  carrier  takes  on  one  of  four  equally  spaced  values 
such  as  0, 71/2,  tc,  and  3n/2,  where  each  value  of  phase  corresponds  to  a  unique  pair  of 
message  bits.  A  QPSK  signal  can  be  defined  as 

=  0<t<T,  1  =  1,23.4.  (19) 

where  Ts  denotes  symbol  duration  (equal  to  twice  the  bit  duration  Tb).  By  using 
trigonometric  identities.  Equation  19  can  be  written  in  the  interval  0  <  f  as: 

W  (0  =  -  l)f ]cos(2)^.()  -  l)f ]sm(2C()  (20) 
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If  the  basis  functions  Oi(0  =  ^2/7^  cosQ-Tiff)  and  02(0  =  ^Jl/Ts  sm(27tfj)  are  defined 
over  the  interval  0<t<T^  then  the  QPSK  signal  can  be  expressed  in  terms  of  the  basis 
signals  as  [Ref.4]: 

V^ir(0  =  {->/^cos[(/-l)|]a>,(f)-V^sin[(i-l)|]O2(0}  /  =  U,3,4.  (21) 

E.  SIGNAL  CONSTELLATION  FOR  QPSK 

The  QPSK  signal  constellation  is  shown  below.  The  constellation  suggests  that  a 
QPSK  signal  can  be  thought  of  as  a  combination  of  two  pairs  of  antipodal  BPSK  signals, 
with  the  two  BPSK  pairs  orthogonal  (in  phase  quadrature)  to  each  other. 
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F.  PROBABILITY  OF  BIT  ERROR  FOR  QPSK 

Since  QPSK  may  be  considered  as  a  combination  of  two  orthogonal  BPSK 

signals,  the  bit  QPSK  error  probability  is  identical  to  the  BPSK  bit  error  probability  given 
by  Equation  18.  Note  however  that  the  QPSK  symbol  error  probability  is  given  [Ref. 4] 
as; 

(22) 

G.  SIMULINK  MODEL  OF  MPSK  DIGITAL  COMMUNICATION  SYSTEM 

In  the  simulation  of  communication  systems,  there  are  two  options  for 
representing  the  modulation/demodulation  process:  passband  and  baseband.  In  passband 
simulations,  the  carrier  signal  is  included  in  the  simulation  model.  By  the  Nyquist 
sampling  theorem,  the  simulation  sampling  frequency  must  be  at  least  twice  the  sum  of 
the  carrier  signal  frequency  and  the  highest  frequency  component  of  the  message  in  order 
to  recover  the  message  correctly.  Since  the  frequency  of  the  carrier  is  usually  much 
greater  than  the  highest  frequency  of  the  input  message  signal,  a  passband  simulation 
requires  a  very  large  number  of  samples  and  excessive  simulation  times.  The  simulation 
time  can  be  dramatically  reduced  by  using  baseband-equivalent  simulation,  also  known 
as  the  low-pass  equivalent  method.  A  baseband  equivalent  model  uses  complex 
envelopes  to  reduce  the  simulation  sampling  rate  and  thus  speed  up  the  simulation.  As 
the  name  suggests,  the  complex  envelope  is  a  complex  function  of  time  that  carries 
information  about  both  the  time-varying  amplitude  and  the  time-varying  phase  of  a 
passband  signal.  Since  the  amplitude  of  an  MPSK  signal  does  not  vary  with  time,  the 


complex  envelope  of  an  MPSK  signal  has  a  constant  magnitude  but  a  time-varying  phase 
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angle  of  the  complex  envelope  for  the  i-th  symbol  is  2jri/M  [Ref.l].  In  this  thesis  the 
baseband  equivalent  model  for  MPSK,  shown  in  Figure  3,  was  developed  and  applied  to 
2PSK  and  4PSK  cases. 

In  order  to  evaluate  the  effect  of  co-channel  interference,  the  model  in  Figure  3 
was  modified  by  adding  another  MPSK  modulator  and  another  message  source. 

The  MPSK  model  has  the  following  main  functional  blocks: 

•  Sampled  Read  from  Workspace 

•  MPSK  modulation  and  MPSK  demodulation  blocks 

•  AWGN  channel 

•  K-step  Delay 

•  Error  Counter 

•  Triggered  Write  to  File. 

In  addition  to  these  blocks,  auxiliary  blocks  such  as  multiplexer,  oscilloscope,  sum,  pulse 
generator,  switch,  counter,  and  so  on  are  also  used. 

The  Sampled  Read  From  Workspace  block  reads  a  row  of  data  for  a  workspace 
variable  (a  matrix)  at  every  data  sampling  point.  The  output  of  this  block  is  a  vector 
whose  length  is  equal  to  the  number  of  columns  of  the  workspace  variable  matrix.  The 
possibility  exists  to  offset  (postpone)  reading  from  workspace  by  a  desired  value  such 
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Figure  3.  Model  of  MPSK  Communication  System 
that  up  to  the  user-defined  offset  time  the  output  of  this  block  is  equal  to  a  user-defined 
initial  value.  When  a  non-zero  time  offset  is  used,  reading  from  workspace  starts  once 
the  simulation  time  exceeds  the  offset.  In  the  case  that  the  last  row  of  the  workspace 
matrix  is  reached  before  completing  the  simulation  (that  is  when  the  number  of  rows  in 
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the  workspace  matrix  is  less  than  the  number  of  simulation  time  instants)  the  first  row  of 
the  workspace  variable  may  be  read  again  if  the  user  has  selected  the  cyclic  read  option 
(otherwise  the  output  is  0). 

The  MPSK  Mod  block  creates  the  complex  envelope  signal  consisting  of  the  real 
and  imaginary  parts  of  the  signal.  This  block  is  a  combination  of  the  MPSK  Map  and  CE 
(Complex  Envelope)  PM  blocks.  The  digital  input  signals  are  in  the  range  [0,M-1],  where 
M  is  total  number  of  different  phases  and  the  phase  shift  for  input  digit  i  is  2ra/M.  The 
MPSK  Mod  output  is  a  unit-magnitude  complex  analog  signal  whose  discrete  values  can 
be  represented  by  M  points  on  the  unit  circle  in  the  complex  plane. 

The  MPSK  Dentod  block  demodulates/demaps  the  complex  input.  This  block  is  a 
combination  of  the  MPSK  Corr  Demod  and  CE  Min/Max  Demap  blocks.  The  MPSK 
Demod  output  is  an  integer  in  the  range  [0,M-1],  where  M  is  the  number  of  phases. 

The  AWGN  (Additive  White  Gaussian  Noise)  Channel  block  adds  white  Gaussian 
noise  with  user-specified  mean  and  variance  to  the  signal  being  transmitted  through  this 
channel  This  block  uses  the  Gaussian  White  Noise  Generator  block  as  the  noise  source. 
The  vector  length  (N)  for  the  initial  seed  entry  determines  the  output  vector  size.  The 
mean  and  the  variance  of  the  noise  are  specified  such  that  the  mean  can  be  either  a  vector 
of  length  equal  to  the  seed  length  (N)  or  a  scalar,  in  which  case  all  the  elements  of  the 
noise  vector  share  the  same  mean  value.  The  covariance  matrix  can  be  one  of  the 
following: 

•  An  N-by-N  positive  semi-definite  matrix  whose  off-diagonal  are  the  correlations 

•  A  length  N  vector,  in  which  case  the  individual  elements  of  the  noise  vector  are 

uncorrelated  but  have  unequal  variances 
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•  A  scalar,  in  which  case  all  the  elements  of  the  noise  vector  are  uncorrelated  but 
share  the  same  variance. 

The  K-Step  Delay  block  delays  its  input  (scalar)  signal  by  K-sampling  time 
intervals.  This  block  function  is  equivalent  to  a  z"*  discrete  time  operator. 

The  Error  Counter  block  detects  the  difference  between  the  port  1  signal  and  the 
port  2  signal  at  the  times  of  the  rising  edge  of  the  port  3  signal  If  the  difference  between 
the  port  1  and  port  2  signals  is  larger  than  a  user-defined  tolerance,  the  error  counter 
increments  its  output  by  1.  The  rising  edge  of  the  signal  at  port  4  resets  the  error  counter 
back  to  0.  The  Error  Counter  block  accepts  scalar  signals  only. 

The  Triggered  Write  to  Workspace  block  writes  the  input  port  1  input  variable  to 
a  workspace  variable  at  the  rising  edge  of  the  input  port  2  trigger  signal.  There  are  six 
entries  to  this  block:  the  workspace  variable  name,  data  type,  the  number  of  trigger  pulses 
between  saved  data,  the  maximum  row  number,  what  to  keep  in  the  case  of  overflow,  and 
the  threshold  value.  The  message  signal  to  the  input  port  1  can  be  either  a  scalar  or  a 
vector  while  the  trigger  signal  to  the  input  port  2  must  be  a  scalar.  The  saved  workspace 
variable  is  a  column  vector  if  port  1  input  is  a  scalar  and  a  matrix  if  the  port  1  input  is  a 
vector.  In  the  later  case  the  first  element  of  the  input  signal  vector  is  saved  in  the  first 
colunm,  the  second  element  of  the  input  signal  vector  is  saved  in  the  second  column,  and 
so  on.  The  block  saves  the  first  record  at  the  first  rising  edge  of  the  trigger  signal  The 
data  can  be  saved  either  as  strings  or  data.  A  pre-defined  number  limits  the  maximum 
row  of  the  output  variable  (this  number  cannot  be  changed  during  simulation). 
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IV.  BIT  ERROR  PROBABILITY  CONVERGENCE 


In  this  chapter  the  convergence  of  the  bit  error  probability  for  BPSK  with 
AWGN  is  investigated.  The  probability  of  bit  error  obtained  by  simulation  converges  to 
the  theoretical  values  as  the  “transmitted”  number  of  symbols  (and  the  simulation  run 
time)  increases.  Convergence  was  tested  for  signal-to-noise  (SNR)  ratios  with  noise 
variances  (the  signal  power  was  kept  constant)  of  0.5,  1,  2,  4,  9,  and  16.  The  signal 
amplitude  is  denoted  A  (A=l),  the  symbol  dmation  is  denoted  T  (T=ls),  the  sampling 
interval  is  denoted  At  (At=T/4=0.25s),  and  the  noise  variance  is  denoted  a  ^ .  With  this 
notation,  the  signal-to-noise  ratio  SNR  in  dB  and  the  theoretical  values  of  the  bit  error 
probability  for  BPSK  are: 

SiVR(^)=101og(-^)  (23) 

20  At 

I  .SNR^ 

Pb=Q(\10  1°  )  (24) 

The  theoretical  values,  calculated  using  Equations  23  and  24,  are  shown  in  Table  H. 
Before  showing  the  convergence  results  for  each  variance/SNR  case,  it  is  useful  to  explain 
the  results  for  one  variance.  For  example.  Table  I  shows  the  convergence  results  for  the 
noise  variance  of  unity.  The  table  has  the  columns  for: 

•  the  number  of  symbols  used  in  the  simulation 

•  the  theoretical  value  of  the  bit  error  probability 

•  the  bit  error  probability  obtained  by  the  simulation  (for  the  corresponding  number 
of  symbols  “transmitted”) ,  and 
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•  the  difference  between  the  theoretical  and  the  simulation  bit  error  probabilities. 
The  table  shows  that  the  simulation  result  converges  to  the  theoretical  result  as  the 
number  of  “transmitted”  symbols  increases. 


Number  of  Symbols 

Theory 

Simulation 

Difference  in  BER 

1000 

0.0230 

0.0250 

8.69% 

10000 

0.0230 

0.0237 

3.04% 

100000 

0.0230 

0.0233 

1.30% 

1000000 

0.0230 

0.0231 

0.43% 

Table  I.  Simulation  Results  for  =  1 

The  estimated  bit  error  probability  as  a  function  of  the  symbol  number  for  each  of  the 
simulations  (from  10^  to  10^  symbols)  is  shown  in  Figures  4  through  7. 
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Figure  6.  BER  Convergence  for  10^  symbols 
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Bit  Error  Probability,  for  1000000 


Figure  7.  BER  convergence  for  10^  symbols 
For  the  case  of  1,000  symbols,  the  simulation  gives  the  estimate  of  the  probability 
of  error  as  0.025,  as  seen  in  Figure  4.  Note  that  each  time  a  simulation  is  repeated  for  the 
same  number  of  symbols  and  the  same  variance  but  with  a  different  seed  value  for  the 
noise  generator,  a  different  result  is  in  principle  obtained.  However,  the  average  of  all 
such  simulation  results  becomes  essentially  a  constant  as  the  number  of  repeated 
simulations  increases.  When  the  number  of  transmitted  symbols  is  increased  to  10,000,  a 


different  estimate  for  the  bit  error  probabifity  is  obtained,  0.0237.  If  the  number  of 
transmitted  symbols  is  increased  to  100,000,  the  bit  error  probability  obtained  is  0.0233 
and  there  is  very  little  change  past  3,000  or  so  symbols.  Lastly,  increasing  the  number  of 
transmitted  symbols  to  1,000,000  yields  the  most  accurate  estimate  for  the  bit  error 
probability  of 0.0231. 

Next,  the  number  of  transmitted  symbols  was  selected  as  10^  and  the  simulation  was  run 
for  the  noise  variances  between  0.5  and  16.  The  results  for  the  estimates  of  the  bit  error 
probability  and  their  difference  relative  to  the  theoretical  bit  error  probabilities  are  shown 
in  Table  H. 


Noise 

Variance 

Theoretical  Pb 

Value 

Simulation  Pb 

Value 

Number  of 

Trials 

Difference 

0.5 

2.339*10'^ 

2.4+10'^ 

1000000 

2.61% 

1.0 

0.0230 

0.0231 

1000000 

0.43% 

2.0 

0.0790 

0.0790 

1000000 

0% 

4.0 

0.1590 

0.1590 

1000000 

0% 

9.0 

0.2520 

0.2529 

1000000 

0.36% 

16.0 

0.3090 

0.3087 

1000000 

0.097% 

Table  H.  Simulation  Results  for  Various  Variances 
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The  convergence  of  the  bit  error  estimates  as  functions  of  the  numbers  of  transmitted 
symbols  is  shown  in  Figures  8  through  12. 
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Figure  8.  BER  Convergence  for  =  0.5 
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Figure  9.  BER  Convergence  for  =  2 
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Figure  12.  BER  Convergence  for  =  16 
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The  bit  error  probability  estimates  obtained  by  simulation  for  various  SNR’s  are  shown  in 
Figure  13  against  the  theoretical  bit  error  probability  for  BPSK  with  AWGN. 


c  1  Qj  SNR  jjg ,  SNR  exp_dB 

Signal  to  Noise  Ratio 


Figure  13.  BER  vs.  SNR 

As  can  be  seen,  the  theoretical  bit  error  probability  and  the  Simulink  results  match  each 
other  closely. 
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V.  SIMULATION  VERIFICATION  FOR  BPSK  AND  QPSK  WITH  AWGN 


In  this  chapter  the  simulation  results  obtained  from  the  (modified)  Simulink  model 
and  a  Matlab  program  will  be  presented  and  discussed.  The  Simulink  model  was  modified 
to  include  co-channel  interference  (discussed  in  the  next  chapter).  In  this  chapter,  only  the 
AWGN  will  be  considered  (the  co-channel  interference  power  was  set  to  0)  for  BPSK  and 
QPSK  modulations.  The  estimates  for  the  bit  error  probability  were  obtained  using  both 
the  Simulink  model  shown  in  Figure  14  and  the  Matlab  program  mpskLm  listed  in  the 
Appendix  M.  The  mpskl.m  program  uses  Communications  Toolbox  functions  rather  than 
Simulink  blocks  to  implement  the  simulations.  Although  a  time-domain  simulation  of  a 
communication  system  is  in  general  possible  using  either  Simulink  block-diagrams  or  the 
equivalent  Matlab  functions,  there  is  a  very  important  difference  between  the  two.  The 
essential  feature  of  Simulink  is  that  simulations  progress  one  step  at  the  time,  that  is  the 
input  data  is  processed  sequentially  (one  at  a  time),  and  are  referred  to  as  “point  data”  or 
“time-flow”  simulations.  On  the  other  hand,  Matlab  simulations  can  process  vectors 
(“blocks”)  of  data  at  each  simulation  step  and  are  thus  referred  to  as  the  “block-data”  or 
“data-flow”  simulations.  Block-data  simulations  are  substantially  faster  than  point-data 
simulations.  However,  block-data  simulations  cannot  be  applied  in  the  time-domain  to 
itnmnlate  feedback-systems  since  such  systems  require  the  responses  to  be  calculated  at 
each  time  instant  before  new  inputs  are  applied. 
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For  both  BPSK  and  QPSK,  only  AWGN  was  considered  and  the  co-channel  interference 
was  made  negligible  by  setting  the  signal-to-interference  ratio  to  +100  dB.  Bit  signal-to- 
noise  (SNR)  ratios  in  the  range  from  -5  dB  to  +12  dB  were  used  in  the  simulations. 

The  simulation  results  were  plotted  using  auxiliary  Mathcad  “scripts”: 

•  plot_sim_noise2.mcd  for  Simulink  BPSK  simulations 

•  plot_mat_noise2.mcd  for  Matlab  BPSK  simulations 

•  plot_sim_noise4.mcd  for  Simulink  QPSK  simulations  and 

•  plot_mat_noise4.mcd  for  Matlab  QPSK  simulations. 

A.  BPSK  BIT  ERROR  PROBABILITY  ESTIMATES 

1.  BPSK  Results  for  the  Simulink  Model 

The  Simulink  model  was  run  from  a  Matlab  script  in  order  to  automate  the  data 
input  and  data  output  operations  for  simulations  run  for  multiple  values  of  signal-to-noise 
ratio  SNR.  The  input  data  file  is  created  using  Matlab  program  prq).m.  Another  Matlab 
program,  sun_psk.m,  reads  the  input  data  file,  starts  the  simulation  (executed  in  a  loop  for 
various  values  of  SNR),  and  saves  the  simulation  output  data.  The  input  data  used  for  the 
simulation  (and  defined  by  running  the  prep.m  Matlab  script)  are: 

Noise  (1)  or  Noise  and  Interference  (2):  1 

Symbol  Duration:  Is 

Oversampling  Factor:  2 

Minimum  Bit  Signal  to  Noise  Ratio:  -5  dB 

Maximum  Bit  Signal  to  Noise  Ratio:  12  dB 

Number  of  values  for  SNR:  10 
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Minimum  number  of  errors  acceptable:  100 
Factor  multiplying  the  error  numbers:  2 
Maximum  size  of  the  random  integer  arrays:  10^ 

File  name  to  save  data:  SIMULINK_NOISE_2 
The  output  data  saved  in  the  ASCII  format  includes: 

•  bit  error  probabilities  (BER) 

•  signal-to-noise  ratios  (SNR) 

•  signal-to-interference  ratios  (SIR)  (set  to  +100  dB  for  the  noise-only  case) 

•  numbers  of  errors  observed 

•  munbers  of  transmitted  symbols. 

The  estimates  (obtained  from  the  Simulink  model)  for  the  BPSK  bit  error  probabihty  are 
shown  in  Figure  15  together  with  the  theoretical  BER  versus  the  SNR . 
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Figure  15.  BPSK  BER  versus  SNR  (Theory  and  Simulink  Model  Estimates) 
2.  BPSK  Results  Obtained  by  Usii^  Matlab  Program  MPSK1.M 
The  same  MPSK  communication  system  that  was  modeled  in  Simulink  using 
Simulink  and  Communications  Toolbox  blocks  was  implemented  in  Matlab  using  Matlab 
and  Communications  Toolbox  functions.  The  same  parameters  were  used  for  the 
simulation: 

Noise  (1)  or  Noise  and  Interference  (2):  1 

Symbol  Duration:  Is 

Oversampling  Factor:  2 

Minimum  Bit  Signal  to  Noise  Ratio:  -5  dB 
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Maximum  Bit  Signal  to  Noise  Ratio:  12  dB 
Number  of  values  for  SNR:  10  or  18 
Minimum  number  of  errors  acceptable:  100 
Factor  multiplying  the  error  numbers:  2 
Maximum  size  of  the  random  integer  arrays:  10^ 

File  name  to  save  data:  MATLAB_NOISE_2 

The  simulation  output  data  were  saved  and  visualized  using  an  auxiliary  Mathcad  script. 
The  theoretical  bit  error  probability  and  the  estimates  obtained  from  the  Matlab  simulation 
are  shown  in  Figure  16  as  functions  of  the  signal-to-noise  ratio  SNR. 
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Figure  16.  BPSK  BER  versus  SNR  (Theory  and  Matlab  Model  Estimates) 

In  this  case,  simulation  results  differ  from  theory.  At  closer  inspection,  one  notices  that 
the  curve  for  the  estimates  is  similar  to  the  theoretical  curve  but  shifted  to  the  right.  Such 
a  shift  corresponds  to  an  SNR  ‘loss”.  Indeed,  if  the  SNR  for  the  theoretical  curve  is 
reduced  by  3  dB,  the  theoretical  curve  and  the  Matlab  estimate  match  each  other,  as 
shown  in  Figure  17.  Therefore,  we  may  conclude  that  the  Simulink  model  gives  very 
accurate  estimates  for  the  bit  error  rate  of  BPSK  with  AWGN  while  the  Matlab  model 
experiences  a  3  dB  SNR  loss  (compensating  for  this  SNR  loss  again  produces  accurate 
estimates).  This  seems  to  indicate  a  fundamental  problem  with  the  Communications 
Toolbox. 
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Figure  17.  BPSK  BER  versus  SNR  (Theory  for  SNR-3dB  and  Matlab  for  SNR) 

B.  QPSK  BIT  ERROR  PROBABILITY  ESTIMATES 
1.  QPSK  Results  for  the  Simulink  Model 


The  Simulink  model  used  for  the  BPSK  is  a  general  model  that  can  be  used  for  any 
MPSK.  Therefore,  the  same  model  was  used  for  QPSK  by  selecting  the  number  of  phases 
M  (an  input  parameter)  as  four.  The  Simulink  model  results  for  the  QPSK  bit  error 
probability  estimates  are  shown  in  Figure  18,  together  with  the  theoretical  bit  error 
probability  and  its  union  bound  given  by  [Ref.4] 
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as  functions  of  the  bit  signal-to-noise  ratio  SNRbit.  We  note  that  the  estimates  differ 
somewhat  from  theory  but  that  this  difference  diminishes  as  the  bit  SNR  increases. 


-10  -5  0  5  10 

bit  sim  ,SNR  yj 

—  —  m  — 

Figure  18.  QPSK  BER  versus  SNR  (Theory  and  Simulink  Model  Estimates) 

2.  QPSK  Results  Obtained  by  Using  Matlab  Program  MPSK1.M 
The  estimates  for  QPSK  BER  obtained  using  Matlab  are  shown  in  Figure  19, 
together  with  the  theoretical  bit  error  probabihty  and  its  union-bound  estimate  given  by 
Equation  25  as  functions  of  the  bit  SNR. 
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Figure  19.  QPSK  BER  versus  SNR  (Theory  and  Matlab  Model  Estimates) 
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VI.  SIMULATION  VERIFICATION  FOR  BPSK  AND  QPSK  WITH 


ADDITIVE  WHITE  GAUSSIAN  NOISE  AND  CO-CHANNEL  INTERFERENCE 

A.  BPSK  BIT  ERROR  PROBABILITY  ESTIMATES 

In  the  preceding  chapter  only  additive  white  Gaussian  noise  (AWGN)  was 
considered  for  BPSK  and  QPSK  simulations.  In  this  chapter  both  AWGN  and  co-channel 
interference  are  considered.  A  second  MPSK  Mod  block  fed  by  an  independent  message 
source  generates  the  co-channel  interference  and  is  added  to  the  sum  of  the  signal  and 
noise.  In  the  Simuluik  model  and  the  Matlab  program  a  range  of  signal-to-noise  ratios 
(SNR)  and  signal-to-interference  ratios  (SIR)  are  used  to  calculate  the  estimates  of  the  bit 
error  probability  as  functions  of  SNR  and  SIR  for  both  BPSK  and  QPSK  The  results  are 
presented  as  families  of  curves,  one  with  SNR  as  the  variable  and  SJR  as  the  parameter 
and  the  other  with  SJR  as  the  variable  as  SNR  as  the  parameter. 

The  parameters  for  the  BPSK  simulations  using  both  the  Simulink  model  and 
Matlab  were: 

Noise  (1)  or  Noise  and  Interference  (2):  2 

The  Number  of  Phases:  2 

Symbol  Duration:  Is 

Oversampling  Factor:  2 

Minimum  Bit  Signal  to  Noise  Ratio:  -5  dB 

Maximiun  Bit  Signal  to  Noise  Ratio:  12  dB 

Minimum  Bit  Signal  to  Interference  Ratio:  -5  dB 
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Maximum  Bit  Signal  to  Interference  Ratio:  12  dB 
Number  of  values  for  SNR:  10 
Number  of  values  for  SJR:  10 
Minimum  number  of  errors  acceptable:  100 
Factor  multiplying  the  number  of  error:  2 
Maximum  si2e  of  the  random  integer  arrays:  10® 

File  name  to  save  data:  NOISE&INTERFERENCE_2 

The  probabilities  of  bit  error  were  saved  as  a  BER  matrix  whose  rows  correspond 
to  constant  values  of  SIR  and  whose  columns  correspond  to  constant  values  of  SNR.  The 
selected  ranges  for  SNR  and  SJR  were  from  -5  dB  to  +12  dB.  For  each  bit  error 
probability  estimate  (each  matrix  entry)  at  least  100  bit  errors  were  observed.  The 
simulations  ran  imtil  the  error  counter  exceeded  the  minimum  specified  number  of  error,  in 
this  case  100.  Upon  exceeding  the  minimum  number  of  required  bit  errors,  the  simulation 
was  restarted  for  another  combination  of  SNR  and  SJR  (to  create  another  BER  matrix 
entry).  High  values  of  SNR  or  SJR  correspond  to  low  noise  and  interference  powers, 
which  in  turn  implies  low  probabihties  of  bit  error.  Therefore,  the  time  required  for  the 
simulation  to  produce  at  least  100  errors  increases  as  the  row  and  column  indices  (SJR 
and  SNR)  increase,  and  the  BER  matrix  entries  in  general  decrease  with  increasing  the 
row  and  column  indices. 

For  QPSK  simulations  the  same  parameter  values  were  used  as  for  BPSK  except 
for  the  number  of  phases  M,  which  are  four  for  QPSK.  The  BPSK  simulation  results  were 
plotted  using  auxiliary  Mathcad  programs  plot_sim_noise&int2.mcd  for  the  Simulink 
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model  and  plot_mat_noise&mt2.mcd  for  Matlab  program  results.  Similarly,  the  QPSK 
simulation  results  were  plotted  using  auxiliary  Mathcad  programs 
plot_sim_noise&int4.mcd  for  Simulink  model  and  plot_mat_noise&int4.mcd  for  the 
Matlab  program. 

1.  BPSK  Results  for  the  Simulink  Model 

The  results  for  the  bit  error  probability  of  BPSK  with  noise  and  co-channel 
interference  obtained  using  the  Simulink  model  are  shown  in  Figure  20  (BER  versus  SNR 
with  SIR  as  parameter)  and  in  Figure  21  (BER  versus  SIR  with  SNR  as  parameter).  The 
increments  for  both  SNR  and  SJR  are  1.889  dB,  starting  from  -5  dB. 
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Figure  20.  BPSK  BER  versus  SNR  with  SJR  as  Parameter  (Simulink  Model  Estimates) 
The  “lowest”  curve  represents  the  theoretical  bit  error  probability  for  BPSK  with  AWGN 
only.  The  increase  of  BER  with  decreasing  SJR  is  evident  from  Figure  20.  The  BPSK 
probability  of  bit  error  as  a  ftmction  of  SJR,  with  SNR  as  a  parameter  starting  at  -5  dB 
and  with  1.889-dB  increments,  is  shown  in  Figure  21. 
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Figure  21.  BPSK  BER  versus  SIR  with  SNR  as  Parameter  (Simulink  Model  Estimates) 
The  “step  function”  represents  the  noise-free  case,  that  is  the  case  of  co-channel 
interference  only.  In  such  a  case  the  probabihty  of  error  is  either  0.5  when  the  co-channel 
interference  power  is  larger  than  the  signal  power  (negative  SJR)  or  0  when  the  signal 
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power  is  larger  than  the  co-channel  interference  power.  (In  order  to  use  the  logarithmic 
scale,  the  value  of  0  for  positive  SJR  has  been  replaced  by  a  small  non-zero  value). 

2,  BPSK  Results  Obtained  by  Usii^  Matlab  Program  MPSK1.M 

The  results  for  the  bit  error  probability  of  BPSK  with  noise  and  co-channel 
interference  obtained  using  Matlab  program  are  shown  in  Figures  22  and  23  (BER  versus 
SNR  with  SIR  as  parameter  starting  at  -5  dB  and  with  1.889  dB  increments)  and  in 
Figure  24  (BER  versus  SJR  with  SNR  as  parameter  starting  at  -5  dB  and  with  1.889  dB 
increments).  Almost  3  dB  difference  can  be  noticed  between  simulation  results  and  the 
theoretical  value.  This  difference  was  explained  in  the  preceding  section. 
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Figure  22.  BPSK  BER  versus  SNR  with  SJR  as  Parameter  (Matlab  Model  Estimates) 


Figure  23  shows  the  BPSK  probability  of  bit  error  as  a  function  of  SJR,  with  SNR  as  a 
parameter.  Again,  the  difference  between  SNR  curves  is  1.889  dB  as  in  SJR  case.  The 
“step  ftmction”  again  represents  the  noise-free  case,  that  is  the  case  of  co-channel 
interference  only. 
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B.  QPSK  BIT  ERROR  PROBABILITY  ESTIMATES 
1.  QPSK  Results  for  the  Simuhiik  Model 

The  results  for  the  bit  error  probability  of  QPSK  with  noise  and  co-channel 
interference  obtained  using  the  Simulink  model  are  shown  in  Figure  24  (BER  versus  SNR 
with  SIR  as  parameter  starting  at  -5  dB  and  with  1.889  dB  increments)  and  in  Figure  25 
(BER  versus  SIR  with  SNR  as  parameter  starting  at  -5  dB  and  with  1.889-dB 
increments). 
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Figure  24.  QPSK  BER  versus  SNR  with  S  JR  as  Parameter  (Simulink  Model  Estimates) 
The  “lowest”  curve  represents  the  theoretical  bit  error  probability  for  QPSK  with  AWGN 
only.  The  increase  of  BER  with  decreasing  SJR  (for  positive  SNR)  is  evident  from  Figure 
24. 
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The  QPSK  probability  of  bit  error  as  a  function  of  SJR,  with  SNR  as  a  parameter 
starting  at  -5  dB  and  with  1.889-dB  increments,  is  shown  in  Figure  25.  The  “step 
function”  again  represents  the  noise-free  case,  that  is  the  case  of  co-channel  interference 
only. 
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Figure  25.  QPSK  BER  versus  SJR  with  SNR  as  Parameter  (Simulink  Model  Estimates) 
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As  can  be  seen  in  Figure  25,  the  curves  for  the  bit  error  probability  tend  to  the  noise-free 
step  function  as  the  SNR  increases.  Also,  the  curves  for  the  negative  SJR  tend  to  0.5, 
mftaning  that  the  bit  error  probability  is  dominated  by  the  interference  for  negative  SJR,  as 
expected. 

2.  QPSK  Results  Obtained  by  Usii^  Matlab  Prt^ram  MPSK1.M 

The  estimates  for  QPSK  BER  obtained  using  Matlab  are  shown  in  Figure  26 
together  with  the  theoretical  bit  error  probability  as  a  function  of  the  bit  SNR  and  with 
SJR  as  parameter  starting  at  -5  dB  and  with  1.889  dB  increments.  Comparing  the  BER 
curves  in  Figures  24  and  26,  we  note  that  they  differ  slightly. 
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Figure  26.  QPSK  BER  versus  SNR  with  SJR  as  Parameter  (Matlab  Model  Estimates) 

The  QPSK  probability  of  bit  error  as  a  function  of  SJR,  with  SNR  as  a  parameter 
starting  at  -5  dB  and  with  1.889-dB  increments,  is  shown  in  Figure  27.  The  “step 
function”  again  represents  the  noise-free  case,  and  the  curves  are  seen  to  tend  to  the  step 
function  as  the  SNR  increases. 
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Figure  27.  QPSK  BER  versus  SJR  with  SNR  as  Parameter  (Matlab  Model  Estimates) 


From  these  figures,  we  see  that  co-channel  interference  acts  like  noise.  When  co-channel 


interference  is  high  (SJR  is  low),  the  bit  error  probability  is  high  and  when  it  is  low  (SJR  is 
high)  the  bit  error  probability  is  low  provided  SNR  is  also  high. 
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From  an  examination  of  Figures  26  and  27,  we  see  that  we  can  maintain  =  10  ^ 
for  SJR<  12  dB  when  SNR>8.22  dB.  Up  to  a  point,  further  increasing  SNR  allows  SJR 
to  further  decrease.  A  similar  result  is  observed  for  BPSK,  as  can  be  seen  from  an 
examination  of  Figures  20  through  23.  Hence,  the  importance  of  AWGN  on  the  effect  of 
co-channel  interference  cannot  be  overlooked. 
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vn.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  objective  of  this  thesis  was  to  model  and  simulate  MPSK  communication 
systems  in  the  presence  of  additive  white  Gaussian  noise  and  co-channel  interference.  The 
models  were  implemented  as  Simulink  block-diagram  models,  using  Simulink  and 
Communications  Toolbox  blocks,  as  well  as  Matlab  programs,  using  Matlab  and 
specialized  Communications  Toolbox  functions.  Baseband  equivalents  of  passband  MPSK 
systems  were  modeled  using  the  complex  envelope  concept.  The  results  for  the  bit  error 
probability  were  obtained  for  BPSK  (2PSK)  and  QPSK  (4PSK)  modulation  types  with 
additive  white  Gaussian  noise  only  and  with  the  sum  of  additive  white  Gaussian  noise  and 
co-channel  interference.  Furthermore,  the  convergence  (with  increasing  numbers  of 
“transmitted”  symbols)  of  bit  error  probability  estimates  obtained  by  simulation  to  the 
theoretical  bit  error  probabilities  has  been  verified  for  several  cases  of  additive  white 
Gaussian  noise  with  different  variances. 

The  simulation  results  For  AWGN  only  are  in  general  very  close  to  the  theoretical 
values  and  the  estimates  of  the  bit  error  rates  converge  to  the  theoretical  values  as  the 
numbers  of  transmitted  symbols  increase.  However,  it  has  been  observed  that  the  Matlab 
model  requires  a  3  dB  SNR  correction  for  BPSK  in  noise  case  only.  With  this  SNR 
correction,  the  simulation  results  for  the  bit  error  probability  match  the  theoretical  results 
and  the  results  for  the  two  implementations  (SIMULINK  and  MATLAB)  match  each 
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other.  The  BPSK  and  QPSK  are  both  sensitive  to  co-channel  interference  almost  to  the 
same  extent. 

As  a  general  rule,  the  sensitivity  of  both  BPSK  and  QPSK  to  co-channel  interference 
can  be  reduced  by  increasing  the  SNR.  Hence,  the  importance  of  AWGN  on  the  effect  of 
co-channel  interference  cannot  be  overlooked.  Of  course,  if  the  SIR  becomes  too  small, 
then  performance  is  poor  regardless  of  the  SNR. 

B.  RECOMMENDATIONS 

This  research  can  be  continued  to  include  forward  error  correction  and  adjacent 
channel  interference  in  the  developed  models.  Since  only  linear  MPSK  systems  were 
considered,  a  future  study  may  address  transmitter  non-linearity.  Theoretical  derivation  of 
bit  error  probability  can  be  obtained  for  MPSK  with  AWGN  and  interference  which  can  be 
either  phase  locked  to  the  signal  or  with  some  given  distribution  (such  as  Uniform  or 
Gaussian)  of  the  phase  difference  between  the  interference  and  the  signal 
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APPENDIX  A.  BATCH_RUN.M  MATLAB  PROGRAM  FOR  CONVERGENCE 


TEST  OF  SIMULATION  USING  BPSK 

%  This  program  runs  multiple  simulations  for  convergence  of  simulink% 
clear 

nxom^symbols  =  input  ('Enter  the  number  of  symbols  [1000]:  '); 
if  isempty (num_syTnbols)  ,  num^symbols  =  1000;  end 


tic 

sim( 'mpsk7_2 ' ,  num_syinbols)  ; 
disp( 'VAR  0.5  Done! ' ) 
toe 

tic 

sim  ( '  mpsk7_l ' ,  num_syinbols )  ; 

disp( 'VAR  1  Done! ' ) 

toe 

tic 

sim{  'mpsk7_2' ,  num_symbols)  ; 

disp ( 'VAR  2  Done ! ' ) 

toe 

tic 

sim(  'mpsk7_4'  ,num_symbols)  ; 

disp ('VAR  4  Done! ') 

toe 

tic 

sim( 'mpsk7_9 ' ,  num_syinbols)  ; 

disp ( 'VAR  9  Done! ' ) 

toe 

tic 

sim(  'mpsk7_16' ,  num_symbols)  ; 

disp ('VAR  16  Done!') 

toe 
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APPENDIX  B.  BER_PLOT.M  MATLAB  PROGRAM  FOR  PLOTTING  THE 
FIGURES  OF  THE  SIMULATION  FOR  EACH  VARIANCES 


%  This  loads  two  files  and  plots  the  result 
clear 

nuin_bits  =  input ('Enter  the  number  of  bits  [1000]:  '); 
if  isempty {num_bits) ,  num_bits  =  1000;  end 
load  d :  Vermiml .  dat 
load  d:\bitniunl.dat 

plot (bitnuml,  errnuml  ./  bitnuml,  'r-') 
grid 

xlabeK'Bit  N\rmber') 

ylabeK'Bit  Error  Probability') 

title ('Bit  Error  Probability  for  variance  1') 

[num_errors  dummy]  =  size(ermuml)  ; 
disp( 'noise  variance=l') 

BER  =  num_errors  /  num_bits 
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APPENDIX  C.  PLOT_SIM_NOISE2  MATHCAD  PROGRAM  FOR  RESULTS  OF 
NOISE  ONLY  CASE  BPSK  USING  SIMULINK  MODEL 


This  program  gives  the  probability  of  bit  error  vs.  signal  to  noise  ratio  using  outputs  obtained 
from  SIMULATION  test  for  noise  only  BPSK. 

BER  :=  READPRN  (’'h:\veys\siimioise2.ber"  ) 

SNR  :=  READPRN  ( ''h:\veys\siimioise2.snr"  ) 


SNR  :=  SNR^ 


m  :=  0 ..  rows  ( SNR  )  -  1 


SNR  jjB  :=-5,-4.9..  13 
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APPENDIX  D.  PLOT_MAT_NOISE2  MATHCAD  PROGRAM  FOR  RESULTS 
OF  NOISE  ONLY  CASE  BPSK  USING  MATLAB  PROGRAM 

This  program  gives  the  probability  of  bit  OTor  vs.  signal  to  noise  ratio  using  outputs  obtained  from 
MATLAB  test  for  noise  only  BPSK. 

BER  :=  READPRN(  "h;\veys\mamoise2.ber"  ) 

SNR  := READPRNC  "h:\veys\mamoise2.snr"  ) 

SNR:=SNR^ 
m:=0..  rows  (SNR)—  1 
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APPENDIX  E.  PLOT_SIM_NOISE4  MATHCAD  PROGRAM  FOR  RESULTS  OF 
NOISE  ONLY  CASE  QPSK  USING  SIMULINK  MODEL 

This  program  gives  the  probability  of  bit  error  vs.  signal  to  noise  ratio  using  outputs  obtained 
from  SIMULATION  test  for  noise  only  QPSK. 

M  :=4 

BER  :=  READPRN(  ••h:\veys\sinmoise4.ber"  ) 

SNR  := READPRN  ( "h:\veys\siiimoise4.siir"  ) 

SNR:=SNR^  m  :=0..  rows  (SNR)  -  1 
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APPENDIX  F.  PLOT_MAT_NOISE4  MATHCAD  PROGRAM  FOR  RESULTS 
OF  NOISE  ONLY  CASE  QPSK  USING  MATLAB  PROGRAM 


This  program  gives  the  probability  of  bit  error  vs.  signal  to  noise  ratio  using  outputs  obtained  frcnn 
MATLAB  test  for  noise  only  QPSK. 

BER:=READPRN[  "h:\veys\matnoise4.ber"  ) 

SNR  :=READPRN[  "h:\veys\matnoise4.snr"  ) 

snr:=snr^ 

m  :=0..  rows(SNR)  -  4 
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APPENDK  G.  PLOT_SIM_NOISE&INT2  MATHCAD  PROGRAM  FOR 
RESULTS  OF  NOISE  AND  INTERFERENCE  CASE  BPSK  USING  SIMULINK 

MODEL 

This  program  gives  the  probability  of  bit  error  vs.  signal  to  noise  ratio  and  probability  of  bit  error 
vs.signat  to  jamming  ratio  using  outputs  obtained  from  SIMULATION  test  for  BPSK. 

BER  := READPRN  ( "h:\veys\smmomt2.ber"  ) 

SNR  :=READPRN(''h:\veys\smmoiiit2.snr'’  ) 

SIR  :=  READPRN  ( "h:\veys\sinmoint2.sjr"  ) 

SNR  :=SNR^  SJR  :=SJR^ 

m  :=0..rows  (SNR)  —  1 

n  :=0..  rows  (SJR)  —  1 

s 

SNRjjg  :=-5..  12 
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B-B 

SNR=INFINnY 

SNR=-5.00dB 

SNR=-3.11dB 

e-e 

SiyR“"1«22  ujd 

SNR=0.67  dB 

SN  R— 2»5u  do 

OT^TT^  A  A  A 

SNR=4.44  oB 

OXTT>  ^ 

SNR=6.33  oB 

. 

SNR=8.22  uB 

SNR=10.11dB 

SNR=12.00dB 
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APPENDIX  H.  PLOT_MAT_NOISE&INT2  MATHCAD  PROGRAM  FOR 
RESULTS  OF  NOISE  AND  INTERFERENCE  CASE  BPSK  USING  MATLAB 


PROGRAM 

This  program  gives  the  probability  of  bit  error  vs.  signal  to  noise  ratio  and  probability  of  bit  error 
vs.signal  to  jamming  ratio  using  outputs  obtained  from  MATLAB  test  for  BPSK. 

BER  :=  READPRN  ( "h:\veys\inatnoint2.ber"  ) 

SNR  .*=READPRN  ( ”li:\veys\inatnoint2.snr”  ) 

SJR  !=  READPRN  ("h:\veys\matnomt2.sjr"  ) 

SNR  :=SNR^  SJR  :=SJR^ 

m  :=0..  rows  (SNR)  —  1 
n  :=0..  rows  (SJR)  —  1 
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,3.872108  10 

-5 


-4 


-2 


0 


M 


-  SJR=INFINITY 

SJR=-5.00dB 
B-B  SJR=-3.11dB 
SJR=-1.22dB 
©-e  SJR=0.67dB 

-  SJR=2.56dB 

-  SJR=4.44dB 

-  S3R=6.33  dB 

-  SJR=8.22dB 

. .  SJR=10.11dB 

- SJR=12.00dB 
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APPENDIX  L  PLOT_SIM_NOISE&INT4  MATHCAD  PROGRAM  FOR 
RESULTS  OF  NOISE  AND  INTERFERENCE  CASE  QPSK  USING  SBMULINK 

MODEL 

This  program  gives  the  probability  of  bit  error  vs.  signal  to  noise  ratio  and  probability  of  bit 
error  vs.signal  to  jamming  ratio  using  outputs  obtained  from  SIMULATION  test  for  QPSK. 

M  :=4 

BER  :=  READPRN  ( ''h:\veys\siimiomt4.ber"  ) 

SNR  :=  READPRN  ("h:\veys\sinmomt4.snr"  ) 

SJR  :=  READPRN  ( "h:\veys\sinmoint4.sjr"  ) 

SNR:=SNR^-3  SJR:=SJR^  BER  ;=.l^  ^  BER 

2  M-  1 

m :=0.. rows  (SNR)  —  1 
n  :s=0..  rows  (SJR )  —  1 
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APPENDIX  J.  PLOT_MAT_NOISE&INT4  MATHCAD  PROGRAM  FOR 
RESULTS  OF  NOISE  AND  INTERFERENCE  CASE  QPSK  USING  MATLAB 

PROGRAM 


This  program  gives  the  probability  of  bit  error  vs.  signal  to  noise  ratio  and  probability  of  bit  error 
vs.signal  to  jamming  ratio  using  outputs  obtained  from  MATLAB  test  for  QPSK. 


M  :=4 

BER  :=  READPRN  ( "h:\veys\matnoint4.ber”  ) 
SNR  :=  READPRN  ( "h;\veys\matnoint4.siir"  ) 

SJR  := READPRN  ( "h:\veys\matnoint4.sjr"  ) 


SNR  :=  SNR^  SJR  :=  SJR^ 


BER  :=BER 


m;=0..  rows  (SNR)  —  1 
n  :=0..  rows  (SJR)  —  1 


HP5 


-  SNR=INF1N1TY 

SNR=-5.00dB 
B-H  SNR=-3.11dB 
SNR=-1.22dB 
e-e  SNR=0.67dB 

-  SNR=2J6dB 

-  SNR=4.44dB 

- SNR=6.33dB 

-  SNR=8.22dB 

-  SNR=10.11dB 

- SNR=12.00dB 
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APPENDIX  K.  PREPAPE.M  MATLAB  PROGRAM  TO  ENTER  THE  SYSTEM 


PARAMETERS  TO  BE  USED  BY  THE  SIMUINK  MODEL 

%%%  This  prepai'es  the  data  file  for  simulink  runs 
clear 

noise_only  =  menu(’  Select: 

’NOISE’,... 

’NOISE  and  INTERFERENCE  ’); 

num_levels  =  input(  Enter  the  number  of  frequencies  M  [2]:  ’); 
if  isempty(num_levels),  numjevels  =  2;  end 

T_S5an  =  input(Enter  the  symbol  duration  T  [1]:  ’); 
if  isempty(T_sym),  T_sym  =  1;  end 

oversampling  =  input(Enter  the  oversampling  factor  [2]:  ’); 
if  isempty(oversampling),  oversampling  =  2;  end 

min_SNR  =  input(Enter  the  MIN  Signal  to  Noise  ratio[-5  dB]:  ’); 
if  isempty(min_SNR),  min_SNR  =  -5;  end 

max_SNR  =  input(Enter  the  MAX  Signal  to  Noise  ratio  [12  dB]:  y, 
if  isempty(max_SNR),  max_SNR  =  12;  end 
if  noise_only  ~=  1 

min_SJR  =  input(Enter  the  MIN  Signal  to  Interference  ratio  [-5  dB]:  ’); 
if  isempty(min_SJR),  min_SJR  =  -5;  end 

max_SJR  =  input(Enter  the  MAX  Signal  to  Interference  ratio[12  dB]:  y, 
if  isempty(max_SJR),  max_SJR  =  12;  end 
end 

if  min_SNR  ==  max_SNR 
num_noise  =  1 ; 
else 

num_noise  =  input(Enter  the  number  of  values  for  SNR  [10]:  0; 
if  isempty(num_noise),  num_noise  =  10;  end 
end 

if  noise_only  ~=  1 
if  min_SJR  ==  max_SJR 
numjam=l; 
else 

num_jam  =  input(Enter  the  number  of  values  for  SIR  [10]:  ’); 
if  isempty(numjam),  numjam  =  10;  end 
end 
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else 

num Jam  =  1; 
end 

min_errors  =  input(!Enter  the  min  number  of  errors  acceptable  [100]:  y, 
if  isempty(min_errors),  min_errors  =  100;  end 

error_factor  =  input(Enter  the  factor  multiplying  the  number  of  errors  [2]:  y 
if  isempty(error_factor),  error_factor  =  2;  end 

initial_num_symbols  =  error_factor*min_errors; 

max_randint  =  input(Enter  the  maximum  size  of  the  random  integer  arrays  [10'^6]: 
if  isempty(max_randint),  max_randint  =  10^6;  end 

file_name  =  input(Enter  the  file  name  to  save  data  [no  ext]:  ’/s’); 

save  sun_data 
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APPENDIX  L.  MPSK.M  MATLAB  PROGRAM  TO  START  SIMULATION 

%%%  This  runs  MPSK  co-channel  interference  with  additive  noise 
clear 

load  sun_data  %This  loads  all  the  input  data 

delta_t  =  T_sym/(2*oversampling) 
seeds  =  randmt(3, 1,1000); 
signaLseed  =  seeds(l); 
noise_seed  =seeds(2); 
interf_seed  =  seeds(3); 

initial_num_symbols  =  error_factor*inin_eiTors; 
tic 

if  num_noise  >  1 

delta_SNR  =  (max_SNR  -  min_SNR)  /  (num_noise  -1); 
else 

delta_SNR  =  0; 
end 

SNR  =  min_SNR  +  [0:num_noise  -  l]*delta_SNR; 

%noise_var_vect  =T_syin/(2*delta_t)  .*  10  (-SNR/10); 

noise_var_vect  =  10  (-SNR/10); 


if  noise_only  ~=  1 
if  numjam  >  1 

delta_SJR  =  (max_SJR  -  min_SJR)  /  (numjam  -1); 
else 

delta_SJR  =  0; 
end 

SJR  =  min_SJR  -i-  [Omumjam  -  l]*delta_SJR; 
interf_gain_vect  =  10  (-SJR/20); 

else 

SJR  =  - 100;  %  There  is  no  Jamming  so  the  SJR  in  dB  is  -infinity 

interf_gain_vect  =  0; 
end 

BER  =  zeros(num_noise, numjam); 
total_symbols  =  initial_num_symbols; 
for  noise_case  =  l:num_noise 
noise_var  =  noise_var_vect(noise_case); 
for  jam_case  =  Lnumjam 
enough_errors  =  0; 
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num_err  =  0; 

num_S)TObols  =  initial_num_symbols; 
rand_int  =  min([num_symboIs  max_randint]); 
total_symbols  =  num_symbols 

while  enough_errors  ~=  1 

interf_gain  =  interf_gain_vect(jam_case); 
rand_int=min([10^6  num_symbols]); 
clear  error_number 

sim(’pskco_bni’,num_symbols) 

[new_errors  err_cols]  =  size(error_number); 
nuni_err  =  num_err  +  new_errors; 
if  num_err  ==  0 

num_symbols  =  num_syinbols*min_eirors; 
randjnt  =  min([num_symbols  max_randint]); 
total_symbols  =  total_symbols  +  num_symbols; 
elseif  (num_err  >  0  &  num_eiT  <  min_eiTors) 

num_symbols  =  (niin_errors  -  num_err)*ceil(total_symbols/num_err); 
total_symbols  =  total_symbols  +  num_symbols; 
else 

enough_eirors  =  1; 
end 

end 

nuniber_of_errors(noise_case,jam_case)  =  num_eiT 
number_of_symbols(noise_case,jani_case)  =  total_symbols 
BER(noise_case,jam_case)  =  num_eiT  /  total_synabols 

end 

eval([’  save  ’file_name  ’.snr  SNR  -asciil); 
eval([’  save  ’file_name  ’.sjr  SJR  -ascii^); 
eval([’save  ’file_name  ’.ber  BER  -ascii’]); 
eval([’  save  ’  file_name  ’.ner  number_of_errors  -ascii]); 
eval([’save  ’file_name  ’.nsy  number_of_symbols  -ascii]) 
end 
toe 
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APPENDIX  M.  MPSK1.M  MATLAB  PROGRAM  FOR  USING  MATLAB 
PROGRAM  OPTION  VICE  SIMULINK  MODEL 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%  Matlab  Version  for  MPSK  %%% 

%%%  Dr.  Jovan  Lebaric  %%% 

%%%  June  14, 1998  %%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%■ 


clear 


%%%%%%%%%%%%%%%%%% 
%%%  User  Input  %%% 
%%%%%%%%%%%%%%%%%% 


j  =  sqrt(-l); 

noise_only  =  menuC  Select: 

’NOISE’,... 

’  NOISE  and  INTERFERENCE  ’);M  =  input(Enter  the  number  of  phases 
(power  of  2)  M  [2]:  ’); 
if  isempty(M),  M  =  2;  end 

T_sym  =  input(Enter  the  symbol  duration  T  [1]:  ’); 
if  isempty(T_sym),  T_sym  =  1;  end 

oversampling  =  input(Enter  the  oversampling  factor  [2]:  ’); 
if  isempty(oversampling),  oversampling  =  2;  end 

min_SNR_bit  =  input(Enter  the  Minimum  BIT  Signal  to  Noise  ratio  in  dB[-5  dB]:  ’); 
if  isempty(min_SNR_bit),  min_SNR_bit  =  -5;  end 

max_SNR_bit  =  input(Enter  the  Maximum  BIT  Signal  to  Noise  ratio  in  dB  [+12  dB]:  ’); 
if  isempty(max_SNR_bit),  max_SNR_bit  =12;  end 

if  noise_only  ~=  1 

inin_SJR_bit  =  input(Enter  the  Minimum  BIT  Signal  to  Interference  ratio  in  dB  [-5 
dB]:  ’); 

if  isempty(min_SJR_bit),  min_SJR_bit  =  -5;  end 
max_SJR_bit  =  input(Enter  the  Maximum  BIT  Signal  to  Interference  ratio  in  dB 
[+12  dB]:’); 

if  isempty(max_SJR_bit),  max_SJR_bit  =12;  end 
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end 


if  max_SNR_bit  ==  min_SNR_bit 
num_noise=  1; 
else 

num_noise  =  input(Enter  the  number  of  values  for  SNR  [18]:  ’); 
if  isempty(num_noise),  num_noise  =  10;  end 
end 

if  noise_only  ~=  1 
if  min_SJR_bit  ==  max_SJR_bit 
numjam=  1; 
else 

numjam  =  input(Enter  the  number  of  values  for  SIR  [10]:  O; 
if  isempty(numjam),  numjam  =  10;  end 
end 
else 

numjam  =  1; 
end 

min_errors  =  input(Enter  the  min  number  of  errors  acceptable  [100]:  y, 
if  isempty(min_errors),  min_err6rs  =  100;  end 

error_factor  =  input(Enter  the  factor  multiplying  the  number  of  errors  [2]:  0; 
if  isempty(error_factor),  error_factor  =  2;  end 

initial_num_symbols  =  error_factor*min_eirors; 

max_randint  =  input(Enter  the  maximum  size  of  the  random  integer  arrays  [lO'^b]:  ); 
if  isempty(max_randint),  max_randint  =  10^6;  end 

file_name  =  input(Enter  the  file  name  to  save  data  [no  ext]:  ’/s’); 


%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%  Parameters  for  MPSK  %%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%%  Symbol  Frequency  %%% 

Fd  =  1/T_sym; 

%%%  Sampling  Frequency  %%% 

Fs  =  oversampling*Fd; 


%%%  Sampling  Interval  %%% 
delta_t  =  l^s; 
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%%%  Seeds  for  Random  Integers  %%% 
seeds  =randint(2, 1,1000); 

signal_seed  =  seeds(l); 
interf_seed  =  seeds(2); 

initial_num_symbols  =  error_factor*niin_errors; 
tic 

%%%  SNR  Array  %%% 
if  num_noise  >  1 

delta_SNR_bit  =  (max_SNR_bit  -  min_SNR_bit)  /  (num_noise  -1); 
else 

delta_SNR_bit  =  0; 
end 

SNR_bit  =  min_SNR_bit  +  [0:num_noise  -  l]*delta_SNR_bit; 
SNR.symbol  =  SNR.bit  +  10*logl0(log(M)/log(2)); 

%%%  Noise  Variance  Array  %%%% 
noise_var_vect  =  10  A  (-SNR_symbol/10); 

%noise_var_vect  =  T_sym/  (2*delta_t)  *  10  (-SNR_symbol/10); 

%%%  SIR  Array  %%% 
if  noise_only  ~=  1 

if  numjam  >  1 

delta_SJR_bit  =  (max_SJR_bit  -  niin_SJR_bit)  /  (numjam  -1); 
else 

delta_SJR_bit  =  0; 
end 

SIR_bit  =  min_SJR_bit  +  [Ornumjam  -  l]*delta_SJR_bit; 
SJR_symbol  =  SJR_bit  +  10*logl0(log(M)/log(2)); 
interf_gain_vect  =  10  (-SJR_symbol/20); 

else 

SJR_symbol  =  100; 

SJR_bit  =  100;  %  There  is  no  Jamming  so  the  SJR  is  large  (infinity) 
end 

BER  =  zeros(num_noise,numJam); 
total_symbols  =  initial_num_symbols; 
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%%%  Noise  and  Interference  Loops  %%% 

for  noise_case  =  1  :num_noise 

noise_var  =  noise_var_vect(noise_case); 

for  jam_case  =  l;num_Jam 

enough_errors  =  0; 
num_eir  =  0; 

num_symbols  =  initial_num_symbols; 
randjnt  =  min([num_symbols  max_randint]); 
total_S3aiibols  =  num_symbols 

while  enough_errors  ~=  1 


randjnt  =  min([10'^6  num_symbols]); 
input_sequence  =  randint(randjnt,  1,[0  M-1],  signal_seed); 


%%%%%%%%%%%%%%%%%%% 

%%%  Modulation  %%%% 

mod_MPSK  =  dmodce(input_sequence,  Fd,  Fs,  ’psk’,  M); 

%%%  Additive  Gaussian  Noise  %%% 
randn(’state’,sum(100*clock)); 

complex_noise  =  randn(randJnt*oversampling,l)*sqrt(noise_var); 
%randn(’state’,sum(200*clock)); 

%complex_noise  =  complex_noise  + 
j*randn(randJnt*oversampling,l)*sqrt(noise_var); 

mod_MPSK  =  mod_MPSK  +  complex_noise; 


%%%  Additive  Interference  %%% 
if  noise_only  ~=  1 

interLgain  =  interf_gain_vect(jain_case); 

mod_MPSK  =  mod_MPSK  +  interf_gain*randint(randJnt*oversampling,  1,[0 
M-1],  interf_seed); 


%%%  Demodulation  %%% 

demod_MPSK  =  ddemodce(mod_MPSK,  Fd,Fs,’psk’,M); 


%%%%%%%%%%%%%% 

%%%  En'ors  %%% 

[symbols_in_error  duinmy_col]  =  find(demod_MPSK  -  input_sequence); 
[new_eiTors  dummy_col]  =  size(symbols_in_eiTor); 
num_err  =  num_err  +  new_errors; 


if  num_eir  ==  0 

num_symbols  =  num_symbols*inin_errors; 
rand_int  =  min([num_symbols  max_randint]); 
total_symbols  =  total_symbols  +  num_symbols; 
elseif  (num_en:  >  0  &  num_err  <  niin_eirors) 

num_symbols  =  (min_errors  -  num_err)*ceil(total_symbols/num_eiT); 
total_symbols  =  total_symbols  +  num_symbols; 
else 

enough_errors  =1; 
end 

end 

number_of_errors(noise_case,jam_case)  =  num_erT 
number_of_symbols(noise_case,jam_case)  =  total_symbols 
BER(noise_case,jam_case)  =  0.5  *  M/(M-1)  *  num_eiT  /  total_symbols 

end 

eval([’  save  ’  file_name  ’.snr  SNR_bit  -ascii’]); 

eval([’  save  ’  file_name  ’.sjr  SJR_bit  -asciil); 

eval([’  save  ’  file_name  ’.ber  BER  -ascii’]); 

eval([’  save  ’file_name  ’.ner  number_of_eiTors  -ascii’]); 

eval([’save  ’file_name  ’.nsy  number_of_symbols  -ascii’]); 
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